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ABSTRACT 

As part of a major program to use isolated Local Group dwarf galaxies as near-field probes of 
cosmology, we have obtained deep images of the dwarf irregular galaxy Leo A with the Advanced 
Camera for Surveys aboard the Hubble Space Telescope. From these images we have constructed a 
color-magnitude diagram (CMD) reaching apparent [absolute] magnitudes of (M475, M 8 i4) > (29.0 
[+4.4], 27.9 [+3.4]), the deepest ever achieved for any irregular galaxy beyond the Magellanic Clouds. 
We derive the star-formation rate (SFR) as a function of time over the entire history of the galaxy. 
We find that over 90% of all the star formation that ever occurred in Leo A happened more recently 
than 8 Gyr ago. The CMD shows only a very small amount of star formation in the first few billion 
years after the Big Bang; a possible burst at the oldest ages cannot be claimed with high confidence. 
The peak SFR occurred wl.5-4 Gyr ago, at a level 5-10 times the current value. Our modelling 
indicates that Leo A has experienced very little metallicity evolution; the mean inferred metallicity is 
consistent with measurements of the present-day gas-phase oxygen abundance. We cannot exclude a 
scenario in which all of the ancient star formation occurred prior to the end of the era of reionization, 
but it seems unlikely that the lack of star formation prior to ss8 Gyr ago was due to early loss or 
exhaustion of the in situ gas reservoir. 

Subject headings: galaxies: dwarf — color-magnitude diagrams — galaxies: evolution 



1. INTRODUCTION 

Dwarf galaxies are the most common class of galaxy 
in the Universe, which gives them significance far be- 
yond their mass. They are simpler in structure and per- 
haps more easily understood than giant galaxies, and 
so their study can shed light on the processes govern- 
ing galaxy evolution. Dwarf galaxies are also important 
to cosmology, because they are similar in mass to pre- 
dictions for the proto-galactic fragments that collapse 
hierarchically to form giant galaxies. One of the most 
interesting testable predictions to come from numerical 
simulations of galaxy formation in a cold dark matter 
(CDM) cosmology is that after the Universe is reionized 
at redshifts z >6, galaxies of total mass <10 8 M Q are pre- 
vented from further accretion of gas by photoionization 
feedback from the cosmic ultraviolet background (e.g., 
iNavarro h Steinmetdfl99l . This may strongly impact 
a dwarf's evolution, but observational evidence is incon- 
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TABLE 1 
Global Properties of Leo A 



Quantity 



Value 



Reference 



Galactic coordinates (£, b) 
distance modulus (m— M)o 
reddening E(B— V) 
Absolute magnitude M 1 ^ 
M(Hi) (10 7 M Q ) 
Total Mass (10 7 M ) 
12+log(0/H) 

Holmberg semi-axes, ag,b[j 
Star-formation rate (M0/yr) 



196.9, +52.4 1 



24.5 ±0.1 


2 


0.021 


3 


-11.7 ±0.2 


4 


1.1 ±0.2 


5 


<20 


5 


7.38' ±0.1 


6 


3f5, 2'.2 


7 


l-2xl0~ 4 


8,9 



Reference s. — ( 1) NED; (2) Dolphin ct al. (200 j); (3 ) 
ISchlegel et all 111991); (4 ) RC 3; (5) lYoung fc Lol iTl996T) : 
(6) Ivan Zee et al.l II2006T); ( 7) IFisher fc Tullyl l!1975fl ~T8) 
IHunter fc Elmegreenl 1120041) ; f9) IJames et alj H2004I) . 

elusive. 

We are fortunate to have a plethora of dwarf galaxies 
available to study in the Local Group; their proximity 
makes it possible to measure their histories in far more 
detail than is possible in even the next-closest groups. 
Because the relative importance of mergers and inter- 
actions as compared to intrinsic properties is not well- 
known, it is important to study dwarfs in as many envi- 
ronments as possible. Therefore we have begun a Hubble 
Space Telescope Advanced Camera for Surveys (ACS) 
program to measure the complete star-formation histo- 
ries (SFH) of six isolated Local Group dwarf galaxies in 
order to search for the imprints of cosmological processes 
on their evolution: the LCID (Local Cosmology from Iso- 
lated Dwarfs) project (Gallart et al., in preparation). By 
measuring several small, isolated galaxies, we will test for 
correlations and patterns in their SFH that may provide 
clues as to the magnitude of the effects of reionization, 
supernova blowout, and other feedback processes. 




Fig. 1. — left: Leo A as seen in a mosaic of SDSS images (courtesy 
D.W. Hogg). The image is 10'5 X 7.'8; North is up and East is to 
the left. Our ACS/WFC field is marked in red and the Holmberg 
radius in cyan, right: 36" X 54" detail of our ACS/WFC field near 
image center. Background galaxies are clearly visible through the 
body of Leo A. The high ratio of bright blue to bright red stars is 
immediately apparent, as is the north-south density gradient. 

Here we present our first results, a new measurement 
of the SFH of the galaxy Le o A (DDO 69). Leo A was 
discovered by iZwickvl (|1942f ) in the course of a search 
for the lowest luminosity galaxies. It is indeed one of 
the le ast luminous gas-rich galaxies known (e.g.. iMateol 
11998ft . The vital statistics of Leo A are given in Table [TJ 
Leo A is a small, blue galaxy conspicuous for its popu- 
lation of youn g, massive stars, with estimated ag es from 
»10 7 -10 8 vr (jDemers et~all Il984t Prolstovi 119961 ). The 
star-formation rate (SFR) has apparently declined since 
that time, although a few small H II regio ns are present. 

Th e first HST observations of Leo A (|Tolstov et all 
Il998f ) revealed a very high fraction of stars aged «1- 
2 Gyr at very low mctallicity, suggesting that only 
~10% of the stellar m a ss wa s contained in older stars. 
iSchulte-Ladbeck et al.1 (|2002f ) obtained slightly deeper 
data in an offset field and measured a lower overall SFR 
and much higher fraaction of ancient stars: they found 
evidence for a median age of w4 Gyr, with the bulk 
of the old stars older than wlO Gyr. Neither of these 
studies reached the necessary depth to convincingly mea- 
sure the fracti o n and age distribution of the oldest stars. 
Dolp hin" et al.l (|2002| ) supplied the first proof of the exis- 
tence of truly ancient stars in Leo A with the discovery 
of RR Lyrae type variables. Our goal here is to provide 
a definitive measurement of how Leo A's SFR has varied 
with time over its entire lifetime. Further details, includ- 
ing spatial patterns within Leo A, comparison to numer- 
ical models, and thorough comparison to other galaxies 
in our sample, will be deferred to future papers, in prepa- 
ration. 

2. DATA ACQUISITION AND PHOTOMETRY 

Leo A was observed with ACS over 16 orbits during 
the time period 26 December 2005 to 8 January 2006. 
Each HST orbit was devoted to a single exposure of 
at least 2400 seconds duration, split to facilitate cosmic 
ray removal. We chose the filter pair F475W (Sloan g') 
and F814W (Broad I) as having the best combination 
of throughput and temperature sensitivity for F-G type 
dwarfs. Our paramount concern in designing the obser- 
vations was to achieve the maximum possible depth of 
photometry, so that our conclusions would be as free as 
possible from difficult-to-characterize systematics. We 
used small dithers between exposures to allow us to cor- 
rect for hot pixels in the detector, and observed when 
the Sun-target angle was over 100°. The total exposure 
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time in (F475W, F814W) was (19,200, 19,520) seconds. 

The Leo A field is shown in Figure [IJ where 
the left panel shows a mosaic of images from 
the Sloan Digital Sky Surve y (fourth data release; 
lAdelman-McCarthy et a l. 2006), centered on the galaxy. 
The footprint of our ACS observation is marked by the 
red square, while the cyan ellipse shows the approxi- 
mate location of the Holmberg radius. It is important 
to remember that a low surface density sheet of red 
giant stars extends ou t to «7f5 (Dol phin et al.1 l2002t 
IVansevicius et al.l l2004V and in the outer reaches of the 
galaxy the ratio of young to old stars must be signif- 
icantly lower than what we find here for the center of 
Leo A. 

The pipeline-processed data were combined using Mul- 
tiDrizzle to eliminate cosmic rays, and ph otometered us- 
ing DOLPHOT, a version of HSTPHOT (lDoh>hinll2000l) 
incorporating the best available ACS calibrations. Ex- 
tended objects and residual hot pixels were rejected 
based on their brightness profiles, leaving 112,000 well- 
measured stars in the color-magnitude diagram (CMD). 
Artificial star tests show that the typical photometric er- 
ror reaches ±0.1 mag for (111475, msw) = (28.7, 27.9); the 
50% completeness limit is at (111475, ni8i4) = (29.0, 27.9). 

The CMD is shown in Figure [2j Where the density 
of stars is high enough to confound the easy estimation 
of the relative number of stars in the various sequences, 
we have plotted contours corresponding to stellar den- 
sity rather than individual stars. The contours cover 
the range from 8-512 stars/decimag 2 , evenly spaced by 
factors of 2. This allows both the overall density dis- 
tribution as well as the fine structure of the principal 
sequences to be shown in a single figure. 

The main features which bear on the measurement 
of the star-formation history, especially at the previ- 
ously ill-constrained earliest times, are the relatively 
bright locations of the peak density of subgiants and 
the main-sequence turnoff; the underdeveloped horizon- 
tal branch and strong, vertically extended red clump; 
and the paucity of upper red giant branch stars. The 
implication of the bright subgiant branch is immedi- 
ately evident by comparison to the location of the 
old, metal-poor isochrone in Fig. [2] 14 Gyr old stars 
would produce a far fainter subgiant branch than ob- 
served. The peak density of subgiants occurs at Mgi4 
w +1.9, more than half a magnitude brighter than 
comparable measure ments in globular clust ers of sim- 
ilar metallicity (e.g.. iRosenberg et al.l l2000l ). but sim- 
ilar to intermedi ate-age clusters in the Small Magel- 
lanic Cloud (e.g., Mighell et al. 1998) and well-matched 
to isochrones with ages ~5 Gyr. This is striking be- 
cause a constant SFR produces a steep and monotonic 
increase in subgiant density wit h age (e.g., Fig. 1 of 
iGallart. Zoccali fc Apar icio 2005]). Because our photom- 
etry reaches the oldest main-sequence turnoff, we can 
unambiguously quantify the galaxy's youth for the first 
time. 

3. DERIVATION OF THE STAR-FORMATION HISTORY 

The quantitative star-formation history of a stellar 
population can be derived from a CMD by comparison of 
the distribution of stars predicted by an ensemble of the- 
oretical models, corrected for the effects of distance, in- 
terstellar reddening, and the completeness and precision 
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Fig. 2.— Our ACS/WFC CMD for Leo A. Individual stars are 
plotted where their density is less than 8 stars/decimag 2 . Contours 
are uniformly spaced by factors of two to show the overall density 
distribution and the fine structure of the stellar sequences. The 
black lines show the main-sequence and lower RGB of isochroncs 
with (Z, age/Gyr) = (0.0001, 14) and (0.001, 5), respectively. 

limit s imposed by the observations (e.g.. lTolstov fc Sahal 
1996). The model predictions are derived from theoret- 
ical isochrones plus a prescription for the intital mass 
function and the influence of unresolved binary stars. 
The adopted isochrone sets have been created in the 
ACS/WFC filter system (L. Girardi, private communi- 
cation), so no transformation of our photometry was re- 
quired. The variation of SFR with time and the evolution 
of stellar metallicity are taken from the combination of 
parameters that results in the best reproduction of the 
CMD. The "best-fit" synthetic CMD is quantified by a 
X 2 -like statistic derived from the definition of the maxi- 
mum likelihoo d for the case of Poisson-distributed data 
(|DolDhinll200l ). 

We adopt a simul ated annealing approach to the min- 
imization problem ( [Kirkpatrick et al.lll983f ). in order to 
find the best CMD without searching all of parameter 
space, while avoiding false local minima. An example 
of the type of local minima encountered in this prob- 
lem is the range of solutions produced by age-metallicity 
degeneracy. Our imple mentation of the proc edure (by 
A.A.C.) is described in ISkillman eTail ((20031 ). O ur re- 
cent tests using various independent algorithms have 
proven that for data of high quality, the details of the 
fitting pr ocedure do not strongly affect the derived SFH 
(see also ISkillman & Gallartl[200l ISkillman et aIll2003L 
and Aparicio & Hidalgo, in preparation). 

The critical parameters (m— M)o and E(B— V) were 
taken from the sources in Table [TJ The so lutions shown 
here w ere obtained usin g isochron e s from iGirardi et alJ 
(2000) and the IMF from lChabrien (|2003h . Extinction in 
the ACS/WFC filter syst em was obtained by interpola- 
tion within the tables of ISirianni et al.l (|2005f ): Agi4 = 



0.039, A 475 = 0.078. 

We solved for the SFH in 9 age bins; the age bins 
are logarithmically spaced to reflect the physical reality 
that the separation between isochrones of different ages 
is a strongly decreasing function of age. Within each age 
bin, a range of metallicities was allowed, from the low- 
est available metallicity t o a value slightly hi gher than 
that corresponding to the Ivan Zee et alJ (|2006l ) metallic- 
ity: 0.0001 < Z < 0.0015. The only constraint on age- 
metallicity relation was to exclude the lowest metallicity 
isochrones from the youngest age bins, as they are seen by 
inspection to be mismatched to the data. To restrict our- 
selves to inferences drawn from the best-modelled phases 
of stellar evolution, we use only the main-sequence and 
subgiant branch in our likelihood calculation. This does 
not significantly change the derived SFH because we in- 
clude >90% of all stars, most of the excluded stars have 
very little age-sensitivity, and their colors are strongly 
model-dependent. 

Our derived history of star formation in Leo A is shown 
in Figure [3] The best-fit SFH is shown for two differ- 
ent age binnings- solid squares and open squares, re- 
spectively. Each point carries its la errorbar in SFR 
and an age errorbar marking the width of the age bin. 
The heavy curve is a cubic spline fit to the results; it 
represents our best estimate of the form of the SFH in 
the central part of Leo A. Because the total number of 
stars is fixed by the data, increasing the SFR in one age 
bin demands a decrease in the adjacent bin(s). Conse- 
quently, the errorbars in Fig.[3]imply that a slightly older, 
broader SFR peak is allowed by the data. The derived 
rates apply to the ACS/WFC field of view, which cov- 
ers 0.6 kpc 2 at a distance of 795 kpc; this is «47% of the 
area within the Holmberg radiu s, and ~29% of the area o f 
Leo A's intermediate-age disk (|Vansevicius et all [20041. 
The redshift scale at top applies to a flat ACDM universe 
with parameters tak en from the WMAP year-3 dataset 
(|Spergel et al.1 120061 ) . The stellar evolution and cosmo- 
logical clocks disagree at early times because of missing 
physics in the stellar models. The most important of 
these is gravitational diffusion settling of light elements, 
which reduces the ages of low-mass stars by 10- 20% 
ijProffitt fc VandenBerd[l99ltlChabover et al.lll992h . 

Both solutions were obtained with fixed values of dis- 
tance modulus, reddening, and binary star properties. 
We tested the sensitivity of the derived solution to rea- 
sonable variations in these parameters, and found that 
for all reasonable parameter choices, the derived SFH has 
the same general form, because of the exceptional photo- 
metric depth and completeness (Cole et al., in prepara- 
tion). Leo A is best described as a late-blooming galaxy, 
in which the great majority of stars are formed after a 
delay of several Gyr. The apparent early peak of star 
formation above a subsequent minimum is not a high- 
confidence feature, as the errorbars indicate. 

The derived metallicity evolution is not plotted in Fig- 
ure [3] because there is almost no trend with age. A 
constant heavy element mass fraction Z = 0.00081q qqq^ 
([Fe/H] as —1.4 ±0.2; ranges represent the rms scatter 
at a given age) is an adequate representation of the re- 
sults, although there is some hint of a trend of increasing 
metallicity with time, from Z = 0.0006 to Z = 0.001. 

4. SUMMARY & DISCUSSION 
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Fig. 3. — The derived SFH of central Leo A. Data points with 
1(7 errorbars show the best fits to the CMD for two age binning 
schemes (solid squares, open squares). The heavy line shows the 
results of a spline fit to the results, our best estimate of the form of 
Leo A's SFH. If isochrone age = lookback time, the redshift derived 
using the third-year WMAP results and flat ACDM cosmology is 
given along the top axis (note that the isochrone age scale and the 
WMAP age of the Universe are incongruent; see text for details). 

As suggested bv lTolstov et all (|1998l ). Leo A is a pre- 
dominantly young galaxy. However, the depth of our 
photometry allows us to much more precisely describe 
its stellar age distribution: within the ACS/WFC field 
of view, 90% of the star formation has occurred more 
recently than 8 Gyr ago. The current star formation 
rate of «10 -4 M Q yr _1 is in good agreement with the 
rates estimated from Ha emission. The large number of 
stars above the 5 Gyr isochrone shows that the SFR in- 
creased dramatically around this time and continued at 
high level for «3 Gyr before declining. However, some 
truly ancient stars are required to reproduce the observed 
CMD and provide the observed RR Lyraes. The 14 Gyr 
isochrone in Fig. [2] traces the lower envelope of subgiant 
stars, confirming that not all of the star formation was 
delayed. 

Integrating o ver the radial profile from 

IVansevicius et all (j2004f ). we find that Leo A only 
astrated ~l-2xl0 6 M Q at ages corresponding to red- 
shifts z > 1. If all of its present-day baryon inventory 
was present at high redshift, then Leo A would have 
been 95% gas and only 5% stars prior to z «2. Al- 
though extreme, the M(H I) /Lb ratio of the young 
Leo A would have been comparable to those of the 
most gas-rich dwarf irregula r syste ms seen today (e.g., 
IWarren. Jerien fc'IC oribalski] |2004f ). Evidently some 
small galaxies can retain or build large gas reservoirs 
without much depletion by star formation. This 
makes Leo A-like objects intriguing candidates for 
proto-galactic fragments that merge into giant galaxies, 
prompting bursts of star formation without leaving 
much stellar residue in their halos. 



Because of uncertainties in assigning absolute ages to 
old stellar populations, and because our time resolution 
decreases with age, we cannot rule out that all of the an- 
cient star formation in Leo A took place prior to reioniza- 
tion. That would require compressing all of the ancient 
star formation in Fig. [3J plus the formation of the un- 
observed Pop III stars that enriched the galaxy to its 
"initial" metallicity [Fe/H] w -1.5, into an -0.5 Gyr 
window. 

Leo A is not the only "late bloomer" in the Local 
Group, but it appears to be an extreme case. An ex- 
ample of a qualitatively simi lar SFH can be found in 
the Leo I (jGallart et all 119991) dwarf spheroidal (dSph) 
galaxy. The two galaxies are of comparable mass, but 
the latter is a gas-free, bound satellite of the Milky Way, 
with plausible triggers for late star formation in tidal in- 
teractions or accretion events. A trigger in the case of 
Leo A is harder to identify. 

The idea that a dwarf galaxy's SFR is tied only to the 
gas infall rate would require an unusual merger history 
for this tiny, isolated system. Leo A would have gained 
virtually all of its gas in the past few Gyr: not enough 
time for it to have gotten to its present remote location 
from the vicinity of M31 or the Milky Way where encoun- 
ters might be expected to be more likely. Given its prob- 
able long history of isolation, it may be that much of the 
gas in Leo A was present from early on, and therefore the 
suppression of early star formation is not attributable to 
the loss or exhaustion of its initial reservoir. It seems 
more plausible that only a small fraction of the H I 
was able to participate in star formation, with the rest 
kept warm in a halo; the UV background or supernova 
feedback are candidate heat sources. Sequestering most 
of the gas in a halo also speeds the chemical evolution 
clock, which is interesting given the prompt enrichment 
and subsequent flat age-metallicity relation in the galaxy. 
The warm halo would have been diffuse and metal-poor, 
resulting in a long cooling timescale and a possible de- 
lay before it could participate in star formation- perhaps 
triggered by a rare infall or interaction event in this iso- 
lated dwarf. 
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